Remark on the core/halo model of 
Bose-Einstein correlations in multiple particle 
production processes * 

K.Zalewski 

M.Smoluchowski Institute of Physics, Jagellonian University, 
Cracow, Reymonta 4, 30 059 Poland, 
e-mail: zalewski@tli.if.uj.edu.pl 
and 

Institute of Nuclear Physics PAN, Cracow 
February 1, 2008 



Abstract 

The core/halo model describes the Bose-Einstein correlations in 
multihadron production taking into account the effects of long-lived 
resonances. The model contains the combinatorial coefficients aj 
which were originally calculated from a recurrence relation. We show 
that aj is the integer closest to the number j\/e. 

PACS 25.75.Gz, 13.65.+i 
Bose-Einstein correlations. 



Bose-Einstein correlations in multiple particle production processes are 
much studied in order to get information about the interaction regions and 
about the hadronization processes. For detailed reviews see e.g. P, j2], 
The starting point is, usually, the factorizeable approximation (see e.g. 
jS]), where the n-particle density matrix for n undistinguishable particles is: 

1 " 

p(pi, • • • , Pn; pi, • • • , Pn) = — II n p(Pjp; Pj-q)' (1) 
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p(p; p') is some single particle density matrix and the summation is over 
all the permutations of the indices of p and p'. The corresponding unsym- 
metrized density matrix is: 



P^(Pl, • • • , Pn; pi, • • • , Pn) = n P(Pp Pj)- (2) 

The quantities usually presented are the n-body correlation functions 

^ / ^ X _ P(Pl,---,Pn;Pl,---,Pn) /o^ 

P^(Pl,---,Pn;Pl,--->Pn) 

As seen from the definitions Cn(p, . . . , p) = n\. There is a number of difficul- 
ties to check this prediction. On the experimental side, a pair of momenta, 



say Pj and p^-, can be reliably measured only if Qij = y—{pi — PjY exceeds 
some Qmin- At present it is difficult to go with Qmin below some (5 — 10) 
MeV (see e.g. |0 and references given there). This is important, because the 
very small Q region is expected to contain narrow peaks due to long-lived 
resonances [7j. There are other difficulties in the small Q region: Coulomb 
corrections, misidentified particles, perhaps effects of coherence. A model 
which assumes that these other factors can be either corrected for or ne- 
glected is the core/halo model^ [Hj, [E!- According to this model for a group 
of identical particles close to each other in momentum space, but not so close 
as not to be resolved experimentally, the single particle density matrix ele- 
ments p{Pjp', p'jq) reach for jg ^ jp a common limit /(p)p(p; p), where p is 
some average momentum of the particles in the group. Of course for jg = jp 
the matrix element in the numerator cancels with the corresponding matrix 
element in the denominator. Thus the correlation function extrapolated from 
the region accessible experimentally to the point p^ — ■ ■ ■ i — Pn is 

C^r''(P,---,P) = E( J )«,/(p)^- (4) 

where aj is the number of permutations of j elements where no element keeps 
its place. 

Let us note the identity 

k 



aj, (5) 



following from the remark that every permutation of k elements can be char- 
acterized by the number j of elements which changed their places and that 



"'^Let us note, however an attempt to include partial coherence into this model 
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the number of choices of these elements is . The formula can be used to 

calculate the coefficient aj when all the coefficients ak with indices k < j are 
known fOl- In this note we derive a simpler formula for the coefficients a„. 
Several derivations of this result can be found in mathematical textbooks. 
Here we use the idea of the proof from ref. 



Let us multiply both sides of(jS))by(— j and sum over k from 
zero to n. On the left-hand side we get 

E 7^ = "!(-!)" (6) 

fc=0 y'' k=0 

and on the right-hand side 



j=0 \ / k=3 



where the second equality follows from the remark that the sum over k yields 
1 for j = n and (1 — l)"^-' = for j < n. Comparing the two sides one finds 



a„, = n'. 



E 



k=0 



The coefficient of n\ tends to when n increases. It is, however, an alter- 
nating series with monotonically decreasing, non-zero terms. For such series 
the sum of the first n elements approximates the limit with an error less than 
the absolute value of the first rejected term. Thus 



e 



<^ (9) 
n + 1 



and an, for n > 0, can be calculated as the integer closest to n\/e. 

The author thanks T. Csorgo for a stimulating exchange of e-mails as 
well as R. Wit and M. Zalewski for discussions on the mathematical aspects 
of the problem. 
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